







(St iC 






















KJOmwtOXlIBWRY 

Pi,cTSaA*i,t^ sc SSI 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 



DATA REDUCTION FOR THE UNSTEADY AERODYNAMICS 
ON A CIRCULATION CONTROL AIRFOIL 

by 

Billy Murel Pickelsimer 
March 1977 

Thesis Advisor: Louis V. Schmidt 

Approved for public release; distribution unlimited. 



1178054 




i 

I 




UNCLASSIFIED 

SeCURITV CLASSIFICATION OF THIS PACE (Whmr^ Dmtm Snfrm¥) 



DUDLEY KNOX LIBRARY 
KAVAL POS'GRADUAft SCKOOl 



REPORT DOCUMENTATION PAGE 



1. nepo^t numien 



READ INSTRUCTIONS 
BEFORE COMPLETING FORM 



2. OOVT ACCESSION NO 



1. PECIPIEMT'S CATALOG NUM8EP 



4. title r«nrf SubfirU; 



Data Reduction for the Unsteady 
Aerodynamics on a Circulation 
Control Airfoil 



5. type OF REPOPT * PERIOD COVERED 

Master's Thesis 
March 1977 



«. PENFONMING one. HE^OMT NUMBEM 



7. AUT HO Rr«; 



t. CONTRACT OR GRANT NUMBERr*) 



Billy Murel Pickelsimer 



*. BENPOMMING OnCANIZATION NAME AnO AOONESS 

Naval Postgraduate School 
Monterey, California 93940 



10. program element. PROJECT, TASK 
AREA 4 WORK UNIT NUMBERS 



II. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 



12. REPORT DATE 

March 1977 



13. NUMBER OF PAGES 



61 



14. monitoring AGENCY NAME A AODNESS^I/ dUi^rwnt from ControlUng OHicm) 



Naval Postgraduate School 
Monterey, California 93940 



IS. SECUAITY CLASS, (ol thia tipett) 



Unclassified 



Is*, declassi^ication/oowngnaoing 

SCHEDULE 



16. Distribution statement (ot c/ii« M^pert) 



Approved for public release, distribution unlimited. 



17. OlSTRrBUTION STATEMENT (ci «/!• •Pmtemet tn Block 20, 1/ dillmrcnt trmn Report) 



IS. supplementary NOTES 



if. KEY WORDS (Continuo on rororoo oipo H nocomoorr onW Idmntiiy by block numkor) 



20. ABSTRACT rConlinuo on rovoroo cidp U nocoocotr 4pW IBomilfr by block m mebe e ) 

Calculating the lift, drag, and pitching moment coefficients 
for an airfoil from the static pressure distribution obtained 
from wind tunnel tests is routine task when steady flow is 
considered, but it is much more complicated when the airfoil is 
operating in an unsteady flow field, similar to that experienced 
by a helicopter rotor blade, produced by an oscillating wind 
tunnel. A data reduction routine capable of condensing the large 



DO , 1473 

(Page 1) 



COITION or I NOV 8S IS OBSOLETE 
S/N 010}.0l4.«e0I I 

1 



UNCLASSIFIED 



SECONITV CLASSI8ICATION Or THIS PAGE (’BHwi Data SutarW) 



UNCLASSIFIED 

it=CU*»lTV CL ASStrtCATlON OF THIS PAGE^H^fi n^tm Cnff ^ 



numbers of data associated with the unsteady investigation, 
as well as a numerical integration algorithm for the unsteady 
aerodynamic coefficients, were developed; however, no unsteady 
data were collected due to hardware failures. The ability 
of the program was demonstrated on previously obtained steady 
and quasi-steady data and sample results were presented. 



DD Form 1473 
^ . 1 Jan 73 
S/N 0102-014-6601 



UNCLASSIFIED 



security CLASS! IflCATION TN»$ ^ACErWh«« Dmtm Enffd) 
2 



Approved for public release, distribution unlimited 



Data Reduction for the Unsteady Aerodynamics 
on a Circulation Control Airfoil 



by 



Billy Murel Pickelsimer 
Lieutenant, United States Navy 
B.S., Oklahoma State University, 1970 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 



from the 

NAVAL POSTGRADUATE SCHOOL 
March 1977 



c.t 






NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



Rear Admiral Isham Linder 
Superintendent 



Jack R, Borsting 
Provost 



This thesis prepared in conjunction with research 
supported by the Naval Air Systems Command 370C. 

Reproduction of all or part of this report is 
authorized . 



Released as a 
Technical Report by: 



ABSTRACT 



Calculating the lift, drag, and pitching moment coeffi- 
cients for an airfoil from the static pressure distribution 
obtained from wind tunnel tests is routine task when steady 
flow is considered, but it^ is much more complicated when the 
airfoil is operating in an unsteady flow field, similar to 
that experienced by a helicopter rotor blade, produced by 
an oscillating wind tunnel. A data reduction routine capable 
of condensing the large numbers of data associated with the 
unsteady investigation, as well as a numerical integration 
algorithm for the unsteady aerodynamic coefficients, were 
developed; however, no unsteady data were collected due to 
hardware failures. The ability of the program was demonstrated 
on previously obtained steady and quasi-steady data and sample 
results were presented. 
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I. INTRODUCTION 



The primary objectives o£ this investigation were: 

1) to report on the data reduction techniques used to ob- 
tain values for both steady and non-steady aerodynamic 
coefficients, resulting from' the experimentally acquired 
static pressure distribution over the model when the tunnel 
was operating at a specified dynamic pressure, with Coanda 
blowing oscillated at various frequencies and amplitudes of 
momentum blowing coefficient; and 2) to present a sample 
of results. Removal of the model from the tunnel for 
modifications in November 1976 precluded any actual process- 
ing of non-steady data consistent with the stated goals . ^ 
However, the data reduction techniques are applicable 
to the pressure distribution obtained with steady Coanda 
blowing and will be demonstrated on such data obtained in 
April 1976. Prior to removal of the model, the tunnel was 
operated without blowing in an effort to generate a set of 
non-steady data to be used in tailoring the data reduction 
procedures to be used downstream in the on-going project. 

In that test, the model experienced self-excited harmonic 
pressures resulting from vortices being shed alternately 



The modifications, to correct for surface delamina- 
tions near the upper jet nozzle, include a steel upper 
surface trailing edge and slot height adjustment to 0.015 
inches. Additional pressure taps were installed on the 
upper surface forward of the Coanda slot, and leaking 
tubes were repaired while the modifications were being 
implemented. 



9 



from the upper and lower surfaces. The manually recorded 
data consisted of the RMS and mean values of the unsteady 
pressure, and the phase angle in respect to a reference for 
each pressure orifice. A more complete description of the 
data acquisition system is given in Section II B. 
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II. EQUIPMENT AND INSTRUMENTATION 



A. AIRFOIL 

The airfoil used in this study was obtained from the 
Lockheed-Calif ornia Company after completion of their 
design feasibility study on circulation control rotors 
(Ref.l). It had an elliptical cross-section with a chord 
of 10.215 inches, thickness ratio of 0.20, and camber of 
3.31 (Figure 1). A slot to facilitate Coanda blowing was 
machined in the upper surface at chord position X/C=0.95, 
with the slot height set at 0.010 inches. Coanda air was 
supplied to the plenum cavity at both ends of the section 
to insure uniform flow along the span. Fifty-three 
pressure taps at mid-span with a chord-wise orientation were 
monitored through two scanivalve transducers for an elec- 
trical representation of the pressure distribution over the 
airfoil. The pressure tap locations are listed in Table I. 

The 2.0 X 2.0 foot unsteady flow wind tunnel in the 
Department of Aeronautics Laboratories at the Naval Post- 
graduate School was used as it is capable of simulating the 
velocity flow field experienced by blade elements of a 
helicopter in forward flight. The model was mounted at an 
angle of attack of - 5 degrees, which closely corresponds to 
the zero-lift angle of attack for the data reported herein. 
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FIGURE 1 

RIRFUIL CRUGS-SEaiDN HITH TAP LUCRT1DN5 



TABLE I 



RIRFOIL PRESSURE THP LDCRTIDNS 



TAP # 

1 

2 
3 
H 
E 
E 
7 
B 
B 

10 
I I 
12 
13 
IH 
IE 
IE 
17 
IB 
IS 
20 
21 
22 
23 
2H 
2E 
2E 
27 



X 

0.000 
0.00S 
0.0EE 
0. 123 
0.21 1 
0.327 
0.E2E 
0.SEH 
1 .HHl 
1 .331 
2.H3E 
2.BEI 
3.3EB 
E.0BE 
B.03E 
7. 132 
7.E27 
B.013 
B.E3E 
3.02B 
S.H30 
3.9EH 
10.070 
10. 1H0 
10. ISI 
10. ISS 
10.21E 



_Y 

0.000 
0.0E0 
0. 13E 
0.220 
0.30E 
0.H00 
0.E1B 
0.710 
0.B7E 
I .004 
1 .0EE 
I . 1B2 
I .31E 
1.3E0 
I .3H3 
1 .230 
I . 1B2 
I .0BB 
0.SE0 
0.0HH 
0.EE2 
0.HH0 
0.3H0 
0.220 
0. 1H7 
0.073 
0 . 000 



TAP » 

2B 

2B 

30 

31 

32 

33 
3H 
3E 
30 
37 
3B 
3S 
H0 
HI 
H2 
H3 
HH 
HE 
HE 
H7 
HB 
HS 
E0 
El 
E2 
E3 



X 

10.200 
10. 1E2 
10. 1 IB 
10.0E1 
3.9H0 
9.B00 
S.E13 
9. 1 IH 
B.EB7 
7.9E0 
7.E70 
7.0S9 
E.033 
E. 1 IE 

H. 0I9 
2.902 
2.H9E 

I . 992 
1 .H7B 
0.9EH 
0.E21 
0.3H9 
0.230 
0.1 IE 
0.03H 
0.010 



Y 

- 0.020 
-0.0H0 
-0. 1E2 
-0.2EE 
-0.372 
-0.H00 
-0.HH3 
-0.E2E 
-0.E90 
-0.E3B 
-0 . EEB 
-0.E7B 
-0.700 
-0.70B 
-0.E9E 
-0.EEE 
-0.EHE 
-0.E10 
-0.EE0 
-0.H90 
-0.3BB 
-0.320 
-0.2E0 
-0. 17B 
-0.09B 
-0.0E2 



X IN INCHES AFT OF LEADINE EDEE 
Y IN INCHES FROM CHBRD LINE 
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B. INSTRUMENTATION 



The data acquisition system evolved from a classic 
manometer board arrangement with readings hand recorded 
(Ref. 2) and reduced to the present system, which employed 
electrical read-outs with an existing Digitec system with 
customized signal conditioning. The digital multiplexing 
circuitry in the Digitec system was interconnected through 
the scanivalve solenoid controllers for automatic scani- 
valve advancing as the output of the digital voltmeter was 
sequentially printed on a paper tape. Pressure sensing was 
accomplished by connecting each of the pressure taps to 
one of the two scanivalves (Figure 2 ) through tubing of 
uniform diameter and length. This was done to insure a 
constant transfer function for all pressure taps, when 
measuring unsteady pressures. Calibration tests using the 
method reported by Johnson (Ref. 3) showed an acceptable 
variation in dynamic gain in the frequency range from 0 to 
300 Hz (Figure 3) , which included the range of practical 
interest. 

For the unsteady aerodynamics analysis (C^=0) , the 
voltage representation of the surface static pressures was 
of the form 

ex=a^+bx(t). (2.1) 

Full analog measuring techniques were employed in measuring 
the mean voltage, RMS voltage, and phase relationship to a 
reference for each pressure tap. The mean voltage was ob- 
tained by processing the signal through a low-pass filter 
circuit (Figure 4) with a two second time constant which 
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corresponds to a corner frequency of 0.08 Hz on a Bode plot. 
Therefore, the influence of the unsteady voltage was essen- 
tially negated and the output readings represented the 
steady components of the signal. 

C. ADVANCED DATA ACQUISITION SYSTEM 

To meet the primary objective mentioned in Section I, 
a more sophisticated data acquisition and reduction system 

2 

will be required and is now being constructed by Englehardt 
at the Naval Postgraduate School in the Department of Aero- 
nautics. The system, based on the Intel MDS-80 micro- 
processor unit, will control the logic of the experiment, 
collect and store data, and handle the entire data reduction 
process. The algorithm for data reduction discussed in 
Section V was proven sufficiently simple and capable of 
performing all required data reduction and numerical integra- 
tion. 



2 

Master of Science thesis by Englehardt with a complete 
description of the entire microprocessor system will be 
published in June 1977. 
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C^=20 jxtd 




Figure 4 

Low-Pass Filter Circuit 
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III. DATA REDUCTION 



A. STEADY DATA PROCESSING 

The aerodynamic coefficients are defined by integral 
equations but must be approximated by numerical integra- 
tions, since data are available only at a finite number 
of pressure tap locations. The approximate steady normal 
force, chord force, and pitching moment coefficients are 



1.0 






0 









Y/C (max) 



^m(TE) 



Y/C (min) 
Y/C (max) 






C3.1) 



•Y/C (min) 

1.0 



(Cp^-Cp^)cVc)d(Vc) + 



0 






Then, with an angle of attack correction for the force 

coefficients and a moment transfer to the quarter-chord 

position, the aerodynamic coefficients are 

C-=C cosa-C sina 
In c 

C,=C sina+C cosct (3.2) 

d n c ^ 

C /"rc-("0.75C . 

m(c/4) m(TE) n 

The pressure data from the April 1976 tests and the 
November 1976 mean pressure data were transferred to a 
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Hewlett-Packard 9830 tape cassette and processed using the 
program listed in Appendix A, which employs the equivalent 
numerical relationship to the equations stated above. 



B. UNSTEADY DATA PROCESSING 

Before the unsteady data were suitable for numerical 
integration, a dynamic gain correction was necessary to 
account for the transfer function discussed in Section II. 

The frequency of the self-induced oscillations was constant 
at 262 Hz when the tunnel airspeed was approximately 100 
ft/s. From Figure 3, the corresponding gain correction, 
jG^I, was approximately 0.5, which implied that the amplitude 
of the surface static pressure was attenuated by a factor 
of one-half at the transducer. All phase angle values were 
relative to tube number 7, which was arbitrarily chosen as 
a clock reference. Phase angle errors resulting from the 
two different scanivalve circuits were faired out in the 
analysis (Figure 5) . 

With the above corrections, the numerical integration 
of the unsteady static pressure distribution was carried 
out in a manner similar to the steady pressure integrations. 
For the unsteady normal force coefficient 



1.0 



(cj) sini> = 



(C sin(J)-C sin4»)d(^/p) 

Pi Pu 



0 




1.0 



f 




(3.3) 



0 
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Then 



<c„^)^ = "^(<Cn^>'sin4.)^C<C„^>\os + ) 



)h 



And 



<|)^ = tan 



-1 ( <l:jj‘)sin+) 

~ n ^ — 

( ( c /} C 0 S 4 .) 



/ 2\'^ 

where <C ) is the RMS value o£ the normal force coefficient 
' n ' 

and (|)^ is the phase angle relative to tube number 7. The 
chord force and pitching moment coefficients were calculated 
similarly. 




Y/C(max) 

' Y 

(C sin(j)-C sin(j))d( /p) 

J '' Pf Pr ^ 'C' 

-Y/C(min) ^ ^ 




Y/C (max) 
-Y/C(min) 



(C^ cos<i)-C 
Pf P 



coS(j))d( /p) 





sin(j)) 




h 

> cos<|>) 



2 






<j>T 




h 

> sincj)) 



> COS(|)) 




Y/C (max) 
-Y/C (min) 



(C sin<i)-C sin(^) (^/p)d(''^/p) +- 

Pf Pj. 



1.0 

^ (C sin,j,-C sin<},)(l-^/.)d(^/p) 

A ^1 
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1 



0 




te 



’ ‘*’t tt? 



2 ^ 

(<Cm > sin4>) 



(<Cm^> coS({)) 



The unsteady aerodynamic coefficients are calculated as in 
Eqs. 3.2, with the aid of elementary trigonometric identities 
to account for the phase angles. The program listed in 
Appendix A has the option for unsteady numerical integration 
incorporated. 

It is intended that the C =0 oscillating pressure 



measurements be repeated when the improved data acquisition 
system becomes available. The tests will be conducted at 
several tunnel dynamic pressures to determine the Reynolds 
number dependence of the data. The previous measurements 



were at approximately R =0.54 x 10^. 
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IV. RESULTS 



A. STEADY DATA RESULTS 

The validity of the concept of circulation control was 
demonstrated by Englar (Refs. 4 and 5) and is confirmed 
here with the steady momentum blowing coefficient results. 
The resulting aerodynamic coefficients are summarized in 
Table II. The ability to produce relatively large in- 
creases in lift coefficient with only a small increase in 
blowing coefficient, shown graphically in Figures 6a, b, 
and c, is one of the desirable characteristics of the 
system. Probably the best measure of the effectiveness is 
the aerodynamic lift amplification ratio, 3C^/3C^, which 
was computed to be approximately 15.3, the low value being 
attributed to the geometric deficiencies of the trailing 
edge slot as alluded to in Section I. However, the results 
are encouraging. 

Close scrutiny of the coefficient of pressure informa- 
tion plotted against X/C (Figures 7-12) reveals a pressure 
distribution over the model such that the center of pressure 
was at the mid-chord position, which is in -agreement with 
the theory. The center of pressure position is also 
verified by the graph of • 




C (3C /3C ) 

= T - nV = 



^ (3C"7”c^J) 



.25 - = .48 



C.P. 



(T3TT) 
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TABLE- II 



AERODYIUMIC COEFFICIENT SUMMARY 
STEADY AERODYNAMIC COEFFICIENTS 
Angle of Attack=-5 degrees U=100 ft/s 



Run No. 


C^ 


=1 




^m(c/4) 


42301 


0 


- 0.0942 


0.0541 


-0.1065 


42502 


0.0855 


1.0352 


0.1085 


- 0.3907 


42503 


0.0640 


0.8146 


0.0678 


-0.3565 


42504 


0.0448 


0.5774 


0.0440 


-0.2776 


42505 . 


0.0272 


0.5459 


0.0327 


- 0.2175 


452O6 


0.0155 


0.1100 


0.0567 


-0.1664 



UNSTEADY AERODYIIAMIC COEFFICIENTS 



Cj^=-0.0765 Cj_(RMS) 

C^= 0.0451 C^(RMS) 

C =-0..1063 C„(RM3) 

m m 



2.1404 


i>=223 


degrees 


0.1629 


II 


degrees 


0.5817 


(|)= 63 


degrees 
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Also, the low drag qualities of the airfoil are favorable 
when full thrust recovery of the momentum blowing coefficient 
is assumed. Pressure coefficient data are listed in 
Appendix B. 

B. UNSTEADY DATA RESULTS 

The unsteady aerodynamics of the elliptically shaped 
airfoil have been given a rather cavalier treatment in 
the past but their contribution to the total forces and 
moments must be realized. The implications of the unsteady 
aerodynamics is not generally well understood; therefore, 
a simple example is presented here to emphasize their 
importance before discussion of the unsteady results. 

Consider an airfoil at zero angle of attack with the 
RMS pressure distribution equal on the upper and lower 
surfaces (Figure 13a) but with unequal relative phase 
angles (Figure 13b) taken from the quarter-chord and the 
three-quarter-chord positions from the actual unsteady data 
(Figure 5). The unsteady normal force coefficient, ) , 

is calculated from the numerical approximation to Eq. 3.3, 
as ; 



) sinij)= 0 . 5 (sin (210°) +sin (230°) -sin (35°) -sin (55°) ) 





-1.33 

-TTTF 




<j> = 223° 
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Figure 13a 



Region 



phi 



1 35 degrees 

2 55 degrees 

3 210 degrees 

4 230 degrees 




Figure I3b 

Sample Pressure Distribution 



Similarly, for the moment coefficient (referenced to the 
trailing edge) : 



Clearly, the force and moment coefficients are of significant 
magnitude even if a uniform distribution (RMS) of pressure 
coefficient were assumed, because of the importance of con- 
sidering relative phase relationships. The dynamic behavior 
of the unsteady aerodynamics of the circulation control 
airfoil, particularly the lift amplification ratio, in the 
unsteady environment caused by cyclic blowing will be in- 
vestigated after the improved data acquisition system is 
operational . 

As discussed in Section I, measurements were taken in 
November 1976, with the airfoil experiencing harmonic 
pressure variations. The unsteady data are listed in 
Appendix C. The experimental data were numerically inte- 
grated, similar to the previous example, for the lift, drag, 
and pitching moment coefficients (Figures 14 and 15) . It 
should be recognized that the phase shift between the upper 




? '2 

(C ) sincj) = -0.6 
' m 

{cj") cos = 0.5((cos(210°)-cos(35°))0.75+(cos(230°)-cos(55°))0.25) 

cos = -0.78 
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and lower surfaces is approximately 180 degrees, the result 
being the relatively large amplitude of the aerodynamic 
coefficients. Highlights of the unsteady results, summarized 
in Table II, were that: 

1. The time varying pressures over the model were 
predominantly harmonic, with a non-dimensional frequency 
(Strouhal no.) of 




2. The RMS value of the normal force coefficient was 
approximately 2.14, which was an order of magnitude greater 
than the RMS drag coefficient. 

3. The test Reynolds no. was approximately 

R = 0.54 X 10^. 
e 

4. The location on the airfoil for the center of 
pressure due to unsteady loadings was approximately 0.491 
chord and the RMS value of pitching moment at this reference 
point was 

(C = 0.1694 

with a phase angle of 90 degrees lag relative to the normal 
force. This term may be considered as an apparent mass 
type effect. 
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V. DATA REDUCTION FOR THE ADVANCED DATA REDUCTION SYSTEM 



The advanced data acquisition system described in 
Section II C will allow a more comprehensive investigation 
of the unsteady aerodynamics of the model, but the data 
reduction techniques must be tailored to handle the large 
numbers of data to be collected. Oscillation of the plenum 
pressure as described by Bauman (Reference 6) induces 
harmonic blowing from the trailing edge slot at a specified 
frequency, W^. It is assumed that the unsteady static 
pressure at each station will have a steady value plus a 
periodic value with frequency superimposed due to the 
oscillating momentum blowing coefficient. 

Since the data are assumed periodic in form, the un- 
steady pressure at any point on the airfoil may be defined 
by a Fourier series (Ref. 7) as follows: 



P(x,t)=A^(x)+J^ [a^(x) cos(nW^t)+B^(x) sin(nW^t)] (5.1) 



where is the frequency of the oscillation. The standard 
formulae for obtaining the Fourier coefficients for the 
infinite series are: 

T 

A_(x) = — [ P(x)dx 

0 2T J 



-T 



Aj^(x) = .j. 



P (x) cos (nx) dx 



(5.2) 



-T 



Bn(3c) j 



T 



-T 



P (x) sin (nx) dx 
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Because the experimental data are in the form of a 
discritization of a truncated periodic signal, the above 
equations are slightly altered for estimating the Fourier 
coefficients. The modified form is: 

V 1 ^ 

(x) = Y ^ 

0 K i-i 1 

. 2 ^ 

A (x) = ^ 2 P. (x) COS (nW^ t) ; n=l,...,n (5.3) 

n ^ i=i ^ ^ 

^ 2 K 

n^^^ = K ^ P^Cx)sin(nw^t) 



where k is the number of samples to be analyzed and n 
is the harmonic number. It is recognized that the physical 
meaning of A^(x) is the steady or mean value of the oscillat 
ing pressure. That is, A^(x) = ^p^(x)) where 



<Po(x)> 



Lim 
T-*- » 




P(x)dx 



(5.4) 



By subtracting out the mean pressure value, the signal becomes 
biased to a pure oscillating pressure with a mean value equal 
to zero. The terms A^(x) and B^(x) are equivalent to the 
cosine and sine components, respectively, discussed in 
Section III B. Figure 16 illustrates the reconstruction of 
a discritized signal utilizing this algorithm programmed on 
the Hewlett-Packard 9830 calculator (program listing in 
Appendix A) . The error was minimized by insuring that the 
number of samples analyzed was an integer multiple of the 
frequency; however, the truncation error was never greater 
than 2%. 
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FIGURE IE 

RECONSTRUCT I UN DF R TRUNCATED SIGNAL 



Having developed an effective method for handling the 
large numbers of data, calculation of the aerodynamic force 
and moment coefficients can be accomplished. Before the 
actual processing of data begins, the sample pressures are 
converted to non-dimensional coefficient form to expedite 
calculations. The unsteady pressure coefficients are 



Cp.(t) = i = 1,2,3,. ..,53 (5.5) 

where P and H are steady pressure values from a pitotstatic 

tube mounted on the floor of the wind tunnel. Then the 

oscillating signal at any pressure tap is; 

C (x,t) = Z A cos(nw,t)+B sin(nw.t) 
p^*-^ n ^I'^n ^1-^ 



Cp(x,t) = 2 T^ sin(nw^t+(|)^) 



(5.6) 



where 



B 



T = '\a ^+B ^ , (j) = n, 

n ^ n n ’ ^n / 



n 



The unsteady Cp(x,t) are related by the phase angles <J> , 
however, the phase angles are meaningless unless they are 
measured relative to a "clock” or reference signal. There- 
fore, it is necessary to designate one pressure tap to be 
used as the reference. Then, the relative phase angle is 



■^rel = '^’i ■ '^’ref 



(5.7) 



For simplicity, shall hereafter be interpreted as the 
relative phase angle. 

With the steady pressure coefficient removed from the 
oscillating signal, separate aerodynamic coefficients can 
be calculated for the steady and unsteady signals. The 
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steady force and moment coefficients are found exactly as 

in Section III A. The unsteady aerodynamic coefficients can 

be defined in a Fourier series form as: 

C^(t) = I sin(kw^t+(j)j^) 

C^(t) = Z sin(kw^t+(j)j^) (5.8) 

k 

C (t) = E C sin (kw^ t+(b, ) 

m^ mj^ '■1 ^k'^ 

The magnitudes and phase angles of the coefficients can be 
determined thusly: 

for the first harmonic 

sin(w^t+(j>) = sin4>cos (w^ t) +C^ cos 4 >sin (w^t) 



where 

C sincf) 
"l 



C COSd) 
"l 



1 K 

0 L i U 

C (-)coS(t., (-) - C (-)cos.j) (-)] d(^) 
I Pi ^l^C*^ p ^C-^ ^U^C-^ '•C^ 



(5.9) 



then 



C 

n 



1 




COS(J>^) 



2 



+ 



(C„^sin*j)2, 



Tan 



sin<t>i 

ni 1 

C cosTT 
"i 1 



Similar calculations . are done to obtain values for the higher 

harmonics (C ,C , ' * * ,C ) as well as the chord force and 
^2 ^3 ^i 

pitching moment coefficients using the higher harmonics of 
pressure coefficient information- as obtained during the 
initial data reduction previously described. 



43 



VI. CONCLUSIONS 



As mentioned in Section I, the two-fold objective of 
this paper was to develop a computer program for unsteady 
aerodynamic analysis and to demonstrate its capabilities 
on experimentally generated unsteady data. The objective 
was partially satisfied; a computer program, well within 
the limitations of mini-calculators such as the Hewlett- 
Packard 9830, was developed which employed a straight- 
forward numerical integration of pressure coefficients to 
obtain values for the aerodynamic coefficients. The pro- 
gram can be modified for unsteady data by properly taking 
into account the cosine and sine components of the pressure 
signal relative to a fixed harmonic (clock) signal. ' Also, 
the extraction of the Fourier coefficients of the truncated, 
time-varying signal appears to be within the capabilities 
of the updated data acquisition system. 
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APPENDIX A 



COMPUTER PROGRAM LISTINGS AND OUTPUT 



1 REM PROGRAM FOR DECOMPOSITION OF A FUNCTION 

2 REM INTO FOURIER COEFFICIENTS 
10 DIM EC2563 

30 W1=(20*PI) 

35 DISP ’’DIVISIONS PER CYCLE?”; 

36 INPUT D 

42 P=2*PI/W1 

43 T=P/D 

50 REM COMPUTE SAMPLE FUNCTION 
60 FOR 1=1 TO 220 
65 X=W1*(I-1)*T 

70 ECi:=SIN(X)+0. 5*SIN(2*X)+0 .3*SIN(3*X) 

80 NEXT I 
100 PRINT 

105 PRINT ’’THE DRIVING FREQUENCY IS 10 HZ” 

106 PRINT D; ’’DIVISIONS PER CYCLE” 

107 PRINT 

110 PRINT ” COEFFICIENTS” 

111 PRINT 

113 PRINT TAB2”COS”TAB15”SIN” 

114 FOR J=1 TO 3 

115 PRINT 

117 PRINT TAB2”A”;J, ”B”;J, ’’PHASE ANGLE” 

119 A1=B1=0 

120 FOR 1=1 TO 220 

130 Bl=(Bl+ECn*SINCJ*Wl*(I-l)*T)) 

140 Al=CAl+ECI3*COS(J*Wl*(I-l)*T)) 

150 NEXT I 

159 P=ATN(B1/A1) 

160 WRITE (15,170) (Al*2/220) , (Bl*2/220) , (P*180/PI) 
170 FORMAT F6 . 3 , 7X ,F6 . 3 , 13X, F5 . 1 

175 NEXT J 
180 END 



THE DRIVING FREQUENCY IS 10 HZ 
10 DIVISIONS PER CYCLE 



COEFFICIENTS 



COS 


SIN 




A 1 


B 1 


PHASE ANGLE 


0.000 


1.000 


90.0 


A 2 


B 2 


PHASE ANGLE 


0.000 


0.500 


90.0 


A 3 


B 3 


PHASE ANGLE 


0.000 


0.300 


90.0 
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1 REM COMPUTE CL, CD, CM- -STEADY AND UNSTEADY 

2 REM FROM CP DATA ON TAPE CASSETTE 
10 COM R1,DC80J ,P[80J ,CC60] 

20 DIM AClOea 
30 PRINT 

40 LOAD DATA 4, A 
60 C9=0 
70 X9=0.95 
80 Y9=0.048 

90 DISP "CP DATA FILE NO. IS"; 

100 INPUT Z 

110 LOAD DATA Z 

111 DC22l= 0.4 

112 DC36J=0.539 

113 Gl=l 

115 PRINT "INPUT OPTION" 

116 PRINT "1 FOR STEADY DATA" 

117 PRINT "2 FOR NON-STEADY DATA" 

118 INPUT G3 

119 GOTO G3 OF 130,212 

120 GOTO 130 
130 PRINT 

140 PRINT "RUN NO. IS"R1 

150 PRINT TAB7"UPPER"TAB24"LOWER" 

160 PRINT TAB2"N"TAB9"CP"TAB19"N"TAB26"CP" 

170 FORMAT F3 . 0 , 2X, F7 . 3 , 5X,F3 . 0 , 2X , F7 . 3 
180 FOR N=1 TO 27 
190 M=N+26 

200 WRITE (15,170)N,C[N] ,M,C[M] 

201 NEXT N 

202 FIXED 3 

203 PRINT "CP1="C9"AT X/C="X9", Y/C="Y9 

204 STANDARD 

205 DISP "CORRECT THE CP VALUES" 

206 STOP 

210 GOTO 270 

212 DISP "INPUT DYNAMIC GAIN"; 

213 INPUT G2 

214 FIXED 3 

216 PRINT "RUN NO. IS"R1 

217 C9 = CQ21] +0.33*(C[223 -CQ21] ) 

218 PRINT "TAP"TAB7"CP(MEAN)"TAB17"CP(RMS)"TAB27"PHI" 
220 FORMAT F3 . 0 , 3X, F7 . 3 , 3X,F7 . 3 , 3X , F4 . 0 

222 FOR 1=1 TO 53 

224 WRITE (15, 220)1, CCIJ,D[I3, PCI] 

226 NEXT I 

227 PRINT "C9="C9 

228 DISP "CORRECT THE CP VALUES" 

230 STOP 

235 GOTO 270 

240 FOR 1=1 TO 53 

242 C[I]=G2*DClJ*SIN(PCIJ/57.296) 

244 NEXT I 

245 C9=(D[21] + 0. 33*(D[2 2] -D[21] )) 

246 Gl=2 

247 GOTO 270 
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250 FOR 1=1 TO 53 

252 C[IJ=G2*D[I]*COS(PClJ/57.296) 

254 NEXT I 

256 C9=(DC21J+0.33*(D[22]-D[21])) 

257 G 1 = 3 

258 GOTO 270 
270 L1=0 

272 L2=0 
274 P1=0 
276 P2=0 

310 FOR 1=1 TO 26 
320 IF 1=21 THEN 490 
330 C1 = -C[I3 
340 C2 = -CCI + 1J 
350 X1=ALU 
360 X2=ACI + 1] 

370 GOSUB 1120 
380 L1=L1+N1 
390 P1=P1+M1 
400 Cl=C[I+2^ 

410 C2=CCI + 26] 

420 Xl=ACI + 27] 

430 X2=A[I + 26] 

440 GOSUB 1120 
450 L2=L2+N1 
460 P2=P2+M1 
470 NEXT I 
480 GOTO 640 
490 C1 = -C[21] 

500 C2=-C9 
510 X1=A[21] 

520 X2=X9 
530 GOSUB 1120 
540 L1=L1+N1 
550 P1=P1+M1 
560 Cl=-C9 
570 C2=-C[22] 

580 X1=X9 
590 X2=A 22 
600 GOSUB 1120 
610 L1=L1+N1 
620 P1=P1+M1 
630 GOTO 400 
640 D1=0 
650 D2=0 
660 P3=0 
670 P4=0 

680 FOR 1=1 TO 13 
690 C1+CCI3 
700 C2 = CCI + 13 
710 Yl=ACl + 533 
720 Y2=ACI + 543 
730 GOSUB 1160 
740 D1=D1+N1 
750 P3=P3+M1 
760 NEXT I 
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770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1205 

1210 

1220 

1225 

1230 

1250 

1260 

1270 

1280 

1290 

1300 



FOR 1=40 TO 52 

C1=C[I3 

C2=C[I+13 

Yl=ACI+533 

Y2=ACI+543 

GOSUB 1160 

D1=D1+N1 

P3=P3+M1 

NEXT I 

FOR 1=14 TO 39 

IF 1=21 THEN 970 

C1=-CCI+13 

C2=C[I3 

Yl=A[I + 543 

Y2=ACI+533 

GOSUB 1160 

D2=D2+N1 

P4=P4+M1 

NEXT I 

GOTO 1200 

Cl=-C9 

C2=C 21 

Y1=Y9 

Y2=A 74 

GOSUB 1160 

D2=D2+N1 

P4=P4+M1 

Cl=-C[22] 

C2=-C9 

Y1=Al75] 

Y2=Y9 

GOSUB 1160 
D2=D2+N1 
P4=P4+M1 
GOTO 950 

Ml=Cl*(X2-Xl)*(l-0. 5*(X2+X13 
M1=M1+(C2-C1)*(X2-X1)*(0.5- (X2/3) - (Xl/63 ) 
N1=0. 5*(C2+C1)*(X2-X13 
RETURN 

N1=0. 5*(C2+C1)*(Y2-Y1) 

M1=0. 5*C1*(Y2-Y1)*(Y2+Y1) 

M1=M1+CC2-C1)*CY2-Y1)*((Y1/6)+(Y2/3)) 

RETURN 

PRINT 

GOTO G1 OF 1210,1500,1600 
PRINT "RUN NO. IS"R1 
PRINT 
FIXED 4 

PRINT "COEFFICIENT BUILDUP IS" 

PRINT TAB2"CNU="L1 
P-RINT TAB2"CNL="L2 
PRINT TAB10"TOTAL CN="(L1+L2) 

PRINT TAB2"CCF="D1 
PRINT TAB2"CCR="D2 
PRINT TAB10"TOTAL CC="(D1+D2) 
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1310 PRINT TAB2”CMU="P1 

1320 PRINT TAB2"CML="P2 • 

1330 PRINT TAB2"CMF=”P3 
1340 PRINT TAB2”CMR="P4 

1350 PRINT TAB10”TOTAL CM=" (P1+P2+P3+P4) "AT T.E.” 
1360 M1=(P1+P2+P3+P4) 

1370 N1=(L1+L2) 

1380 C1=(D1+D2) 

1390 M5=M1- CO. 75*N1) 

1400 DISP "INPUT ANGLE OF ATTACK (DECS.)"; 

1410 INPUT A2 

1415 A2=0.01745*A2 

1416 A9=A2 

1420 N5=N1*C0S(A2)-C1*SIN(A2) 

1430 C5=C1*C0S(A2)+N1*SIN(A2) 

1440 PRINT 
1445 C2=A2*57.296 

1450 PRINT TAB 2"CL (MEAN) ="N5 ;TAB21"CD(MEAN) ="C5 ; 
1460 PRINT TAB2"ALPHA ="C2"DEGS" 

1470 STANDARD 

1480 GOTO G3 OF 90,240,250 
1490 END 
1500 N3=L1+L2 
1502 C3=D1+D2 

1504 M3+P1+P2+P3+P4 

1505 FIXED 4 

1506 PRINT "CNCRMS)*SIN(PHI)="N3 
1508 PRINT "CCCRMS)*SIN(PHI)="C3 

1510 PRINT "CM(RiMS)*SIN(PHI)="M3 

1511 PRINT 
1595 GOTO 250 
1600 N4=L1+L2 
1602 C4=D1+D2 

1604 M4=P1+P2+P3+P4 

1605 FIXED 4 

1606 PRINT "CN(RMS)*COS(PHI)="N4 
1608 PRINT "CC(RMS)*COS(PHI)="C4 

1610 PRINT "CM(RMS)*COS(PHI)="M4 

1611 PRINT 

1612 Nl=SQR(N3+2+N4i2) 

1613 P1=ATN(N3/N4) 

1614 X=N3 

1615 Y=N4 

1616 F=P1 

1617 GOSUB 1900 

1618 P1=F 

1620 PRINT "CN(RMS)="N1,"PHI="P1*57.296 
1630 Cl=SQR(C3i2+C4i2) 

1640 P2=ATN(C3/C4) 

1641 X=C3 

1642 Y=C4 

1643 F=P2 

1644 GOSUB 1900 

1645 P2=F 

1650 PRINT "CC(RMS)="C1,"PHI="P2*57. 296 
1660 Ml=SQR(M3+2+M4+2) 
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1670 P3=ATN(M3/M4) 

1671 X=M3 

1672 Y=M4 

1673 F=P3 

1674 GOSUB 1900 

1675 P3=F 

1680 PRINT "CM(RMS) = '’M1,"PHI=''P3*57.296 
1690 A1=N1*C0S(A9)*SIN(P1)-C1*SIN(A9)*SIN(P2) 

1700 B1=N1*C0S(A9)*C0S(P1)-C1*SINCA9)*C0S(P2) 

1710 L9=SQR(Al+2+Blf 2) 

1720 L8=ATN(A1/B1) 

1721 X=A1 

1722 Y=B1 

1723 F=L8 

1724 GOSUB 1900 

1725 L8=F 

1730 A2=N1*SIN(A9)*SIN(P1)+C1*C0S(A9)+SIN(P2) 

1740 B2=N1*C0S(P1)*SIN(A9)+C1*C0S(A9)*C0S(P2) 

1750 D9=SQR(A2+2+B2+2) 

1760 D8=ATN(A2/B2) 

1761 X=A2 

1762 Y=B2 

1763 F=D8 

1764 GOSUB 1900 

1765 D8=F 

1770 A3=Ml*SIN(P3)-0. 75*N1*SIN(P1) 

1780 B3=M1*COS(P3)-0. 75*N1*C0S(P1) 

1790 M9=SQR(A3+2+B3+2) 

1800 M8=ATN(A3/B3) 

1801 X=A3 

1802 Y=B3 

1803 F=M8 

1804 GOSUB 1900 

1805 M8=F 
1810 PRINT 

1820 PRINT TAB10”AERODYNAMIC COEFFICIENT SUMMARY” 
1830 PRINT 

1840 PRINT ”CL=”N5" +”L9”SIN(W1*T+”L8*57. 296”)” 

1850 PRINT 

1860 PRINT ”CD=”C5”+”D9”SIN(W1*T+”D8*57. 296”)” 

1870 PRINT 

1880 PRINT ”CM(C/4)=”M5”+”M9”SIN(W1*T+”M8*57. 296”)” 
1890 GOTO 90 

1900 IF X<0 AND Y<0 THEN 1920 

1910 IF X>0 AND Y<0 THEN 1920 

1915 GOTO 1930 

1920 F=F+PI 

1930 RETURN 

1940 STOP 
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STEADY DATA OUTPUT 



RUN NO. IS 123 

COEFFICIENT BUILDUP IS 
CNU= 0.3480 
CNL=-0. 4281 

TOTAL CN=-0.0802 
CCF=-0.0339 
CCR= 0.0721 

TOTAL CC= 0.0382 
CMU= 0.1092 
CML»-0.2840 
CMF- 0.0017 
CMR= 0.0038 

TOTAL CM=- 0.1693 AT T.E. 

CL(MEAN)=-0.0765 CD(MEAN)= 0.0451 (CM) C/4 (MEAN) =-0 . 109 2 

ALPHA =-4.9991 DECS 



UNSTEADY DATA OUTPUT 



CN(RMS)*SIN (PHI) =-1.5385 
CC(RMS)*SIN (PHI) =-0.0312 
CM(RMS)*SIN (PHI) =-0.6336 

CN (RMS) *COS (PHI) =-1.4968 
CC(RMS)*COS(PHI)=-0.0042 
CM(RMS)*COS (PHI) =-0. 8581 



CN(RMS)= 2.1465 
CC(RMS)= 0.0315 
CM(RMS)= 1.0667 



PHI=225.7868 
PHI=262. 2760 
PHI=216.4451 



AERODYNAMIC COEFFICIENT SUMMARY 



CL=-0.0765 


+ 2.1405 


SIN(W1*T+ 


225.8305 


) 


CD= 0.0451 


+ 0.1629 


SIN(W1*T+ 


39.2111 


) 



CM(C/4)=-0. 1092 + 0.5836 SIN(W1*T+ 63.0435 ) 
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RPPENDIX B 

PRESSURE CQEFFICIENT DPTP WITH STEADY <;^ 

RUN ND. H2321 Oc=0.00 

RNELE DF RTTflCK=*-E DEE U=100 FT/SEC 



PRESSURE PRESSURE 

TAP » CQEFFICIENT TAP » CDEFFICIENT 



1 


0.230 


2B 


-0.0B0 


2 


0.E37 


23 


-0. 1 IE 


3 


0.32B 


30 


-0.0ES 


H 


0.BB7 


31 


-0. 12B 


E 


0.B1B 


32 


-0. 107 


E 


0.70B 


33 


-0. 133 


7 


0.HH3 


3H 


-0. 171 


B 


0. lEE 


3E 


-0. IBB 


3 


0.00B 


3B 


-0.201 


10 


-0. 137 


37 


-0.2B1 


1 1 


-0.2E3 


3B 


-0 . 2B2 


12 


-0.3E3 


33 


-0 . 2BH 


13 


“0.H3H 


H0 


-0.333 


IH 


-0.B1 1 


HI 


-0.3HB 


IE 


-0.EEE 


H2 


-0.331 


IE 


-0.E0H 


H3 


-0.E03 


17 


-0.E32 


HH 


-0.E03 


IB 


-0.E31 


HE 


-0.EE3 


13 


-0.HBB 


HE 


-0.E72 


20 


-0.HBB 


H7 


-0.B7E 


21 


-0.330 


HB 


-1 .32E 


22 


-0.2HH 


H3 


-1 .333 


23 


-0.230 


E0 


-0.322 


2H 


-0.07E 


El 


-1 .H02 


2E 


-0. IHH 


E2 


-0.3H0 


2B 


-0. IBH 


E3 


-0. 13E 


27 


-0. 1 IB 
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PRESSURE COEFFICIENT DHTH WITH STEADY 

RUN NO. H2302 G^«0.0BES 

HNBLE DF flTTRCK=-E DEE 0=100 FT/SEC 



PRESSURE PRESSURE 

TRP # CDEFFICIENT THP # CDEFFICIENT 



1 


1 .000 


2B 


0. 1 IE 


2 


0.37H 


23 


-1 . 171 


3 


0.E33 


30 


-1.072 


H 


0.H71 


31 


-0.50H 


E 


0.233 


32 


-0.373 


E 


0.033 


33 


-0.253 


7 


-0. IBE 


3H 


0.013 


B 


-0.HHB 


35 


0.05E 


B 


-0.E3B 


3B 


0.033 


10 


-0.B1B 


37 


0.0E1 


I 1 


-0.35B 


3B 


0.023 


12 


-1 .07H 


33 


-0. 105 


13 


-1 .20S 


H0 


-0.100 


IH 


-1 .27H 


HI 


-0.03E 


IS 


-1 .3H3 


H2 


-0. 172 


IE 


-1 .3S5 


H3 


-0.233 


17 


-1.372 


HH 


-0.221 


IB 


-1 .317 


H5 


-0 . 253 


IS 


-1.235 


HE 


-0.215 


20 


! 

LD 


H7 


-0.273 


21 


-3.371 


HB 


-0.3B3 


22 


-5.552 


H3 


-0. 177 


23 


-2.B1H 


50 


-0.020 


2H 


-0.HEB 


51 


0. 1 13 


2S 


-! . I5H 


52 


0.H53 


2B 


-2.7IB 


53 


0.BI 1 


27 


-2 . 3HB 
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PRESSURE CDEFFICIENT PATH HITH STEADY 

RUN ND. M2303 Cu=0.0EH0 

FINELE nr HTTHCK=-E DEB U=100 ET/SEC 



PRESSURE PRESSURE 

TRP # CDEFFICIENT TRP # CDEFFICIENT 



I 


I .000 


2B 


0.339 


2 


0.3B7 


2S 


0,203 


3 


0.7B2 


30 


-0.B9S 


H 


0.5HI 


31 


-0.279 


S 


0.323 


32 


-0.22B 


E 


0. ISE 


33 


-0. I5I 


7 


-0.051 


3H 


0.005 


B 


-0.3B7 


35 


0.0IH 


3 


^0.HBH 


3E 


0.023 


10 


-0.E95 


37 


-0.023 


I I 


-0.BI5 


3B 


-0.031 


12 


-0 . S3 1 


39 


-0.097 


13 


-I .0B5 


H0 


-0. 1 IB 


IH 


-I . IE0 


HI 


-0. IH0 


IS 


-I .235 


H2 


-0.22B 


IE 


-I .23E 


H3 


-0.2EB 


17 


-1 . 155 


HH 


-0.31 I 


IB 


-I .210 


H5 


-0.295 


IS 


-I . 125 


HE 


-0.313 


20 


-I. 051 


H7 


-0.3B7 


2! 


-2.5B3 


HB 


-0.5BH 


22 


-H. 1 m 


H9 


-0.BB7 


23 


-I .BB0 


50 


-0.231 


2H 


-0.2HB 


5! 


0.037 


25 


-0.71B 


52 


0.1 17 


2S 


-2.0HE 


53 


0.709 


27 


-1 . 151 
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ii 




PRESSURE CDEFFICIENT DRTR MITH 5TEH0Y U 



RUN NU. H230H g.«0.0HHS 

RNELE UP RTTFICK=-E DEB U=I00 FT/SEC 



PRESSURE PRESSURE 

TRP # COEFFICIENT TAP # COEFFICIENT 



I 


0.B72 


2B 


0.237 


2 


I .00S 


29 


0. 192 


3 


0:B03 


30 


“0.23B 


4 


0.EBS 


31 


“0.0B7 


S 


0.S3S 


32 


“0. I3E 


S 


0.312 


33 


“0. IB4 


7 


0.030 


34 


“0.0IB 


B 


“0. 197 


35 


“0.03E 


S 


“0.3B2 


3B 


“0.054 


(0 


“0.513 


37 


“0.07E 


1 1 


“0.BI9 


3B 


“0. 105 


12 


“0.775 


39 


“0. I 19 


13 


“0.92B 


40 


“0. 143 


IH 


“I .010 


41 


“0. IB7 


IS 


“I .092 


42 


“0.294 


IB 


“I .02B 


43 


“0.3E2 


17 


“I .054 


44 


“0.3B7 


IB 


“I .0IB 


45 


“0.3E7 


13 


“0.BB7 


4E 


“0.3E7 


20 


“0.903 


47 


“0.534 


21 


“1.949 


4B 


“0.777 


22 


“2.994 


49 


“0.EI3 


23 


“1.034 


50 


“0.297 


2H 


“0.045 


51 


“0.379 


2S 


“0.229 


52 


0.030 


2S 


“0.7B2 


S3 


0.579 


27 


“0.00B 
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PRESSURE CDEFF 


I Cl ENT 


DRTR W 


ITH STEADY Qix 




RUN ND. 42305 
ANGLE OF RTTRCK=“5 


DEE 




C>=0.0272 
U=100 FT/SEC 


TRP 


PRESSURE 

# CDEFFKIENT 




TRP 


PRESSURE 

# COEFFICIENT 




1 


0.723 




2B 


0. 1 17 


2 


0.3B0 




23 


0. 154 


3 


0.3S3 




30 


0.050 


H 


0.aiE 




31 


0.045 


S 


0.E02 




32 


-0.032 


G 


0.423 




33 


-0.0B2 


7 


0.2Ea 




34 


-0.042 


a 


-0.0a7 




35 


-0 . 0E5 


3 


-0.251 




3G 


-0.0B7 


10 


-0.4ia 




37 


-0. I2B 


1 1 


-0.503 




3B 


-0. IBE 


12 


-0.E04 




33 


-0.224 


13 


-0.75E 




40 


-0.224 


IH 


-0 . B37 




41 


-0.224 


15 


-0.3ia 




42 


-0.230 


)G 


-0.a33 




43 


-a.33E 


17 


-0 . BEE 




44 


-0 . 432 


la 


-0.BS1 




45 


-0.451 


13 


-0.E31 




4G 


-0.44B 


20 


-0.B01 




47 


-0.E03 


21 


-1 .370 




4B 


-0.aaa 


22 


-! .340 




43 


-0.B42 


23 


-0.53E 




50 


-0.743 


2H 


0.057 




51 


-0.533 


25 


-0.04E 




52 


-0.2E2 


2G 


-0.030 




53 


0. 1B3 


27 


0.0EE 









56 





PRESSURE COEFF 


ICIENT 


DBTB W 


ITH STEBDY 




RUN ND. H230E 
BNBLE QF BTTRCK=-5 


DEE 




C>i»0.0133 
U=100 FT/SEC 


TBP 


PRESSURE 

# CDEFFICIENT 




TBP 


PRESSURE 
* CDEFFICIENT 




! 


0.522 




2B 


0.0BE 


2 


0.90B 




29 


0.0B5 


3 


0.9EH 




30 


0.075 


H 


0.B2S 




31 


0.0B3 


5 


0.7H3 




32 


0.039 


E 


0.590 




33 


-0.039 


7 


0.31B 




3H 


-0 . 0B0 


B 


0.0HE 




35 


-0. 150 


3 


-0. 13E 




3B 


-0.220 


10 


-0.31 1 




37 


-0.279 


1 1 


-0.H07 




3B 


-0.27E 


12 


-0.HEB 




39 


-0.251 


13 


-0.E0H 




H0 


-0.291 


IH 


-0.E97 




HI 


-0.331 


15 


-0.7B9 




H2 


-0.H07 


IE 


-0.7E2 




H3 


-0.515 


17 


-0.B50 




HH 


-0.579 


IE 


-0.E99 




H5 


-0.5E3 


19 


-0.E32 




HE 


-0.B7H 


20 


-0.EH1 




H7 


-0.772 


21 


-0.975 




HB 


-1 .003 


22 


-1 .309 




H9 


-1 . 123 


23 


-0.30E 




50 


-0.H5H 


2H 


0.00E 




51 


-1 .075 


25 


0.0E7 




52 


-0 . B3B 


2E 


0.097 




53 


-0.290 


27 


0.0E7 
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RPPENDIX C 

UN5TEHDY PRESSURE DPTH 



PNBLE DF RTTPCK=-S DEB U=I00 FT/5EC 

MERN PRE55URE UN5TERDY PRESSURE PHASE 

TRP # COEFFICIENT COEFFICIENT RNELE 



1 


0.270 


0.H0H 


1E0 


2 


0.725 


0. IBE 


1E0 


3 


0.90B 


0.079 


290 


H 


0.B12 


0. IBB 


330 


S 


0.750 


0.203 


3H5 


E 


0.B2H 


0.210 


353 


7 


0.3B3 


0.225 


0 


B 


0.1B1 


0.3BH 


17 


S 


-0.0HB 


0.H00 


23 


10 


“0.1BB 


0.5H5 


32 


1 1 


-0.273 


0.513 


3H 


12 


-0.352 


0.BB1 


37 


13 


-0.HB0 


0.730 


HE 


IH 


-0.553 


0.737 


H9 


IS 


-0.E2E 


0.7H3 


51 


IE 


-0.E21 


0.7B0 


53 


17 


-0.595 


0.E13 


55 


IB 


-0.597 


0.EE1 


E0 


19 


-0.H9B 


0.500 


EB 


20 


"0.H21 


0.5B2 


75 


21 


-0.3E7 


0.HB1 


BB 


22 


-0.313 


0.H00 


103 


23 


-0 . 295 


0.33B 


1 1 1 


2H 


-0.2B3 


0.2E5 


109 


25 


-0.275 


0. IBS 


92 


2E 


-0.33B 


0.1BB 


B9 


27 


-0.H25 


0. laa 


E5 
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UNSTERDY PRESSURE DRTfl 





(CDNT 


INUED) 




RNELE DF flTTflCK=-S DEE 


U=100 FT/SEC 






MEHN PRESSURE 


UNSTE0DY PRESSURE 


PH0SE 


TRP # 


CDEFFICIENT 


CDEFFICIENT 


RNELE 




20 


“0.H1H 


0. 125 


3H0 


20 


-0.205 


0.09H 


295 


30 


-0.230 


0. lEl 


290 


31 


-0.257 


0.223 


290 


32 


-0. 101 


0.31H 


290 


33 


-0. 173 


0.37B 


203 


3H 


-0. 170 


0.H21 


270 


3S 


-0. 192 


0.H70 


255 


3B 


-0.205 


0.539 


2H0 


37 


-0.220 


0 . E00 


237 


30 


-0.2H3 


0.500 


230 


39 


-0.2E0 


0.000 


231 


H0 


-0.20H 


0.B0H 


227 


HI 


-0.300 


0.000 


22H 


H2 


-0.3E0 


0.013 


217 


H3 


-0.500 


0.013 


215 


HH 


-0.550 


0. 000 


210 


HE 


-0.593 


0.530 


205 


HE 


-0.BH0 


0.513 


195 


H7 


-0.053 


0.H03 


190 


H0 


-1 .250 


0.H25 


103 


H9 


-1.300 


0.H05 


175 


S0 


-0.903 


0.H05 


171 


SI 


-! .H53 


0.H50 


152 


52 


-1 . 125 


0.500 


157 


53 


-0.335 


0.H35 


1H2 
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